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ABSTRACT 


This is a proposa] for an experiment to measure the relative 
stopping powers of Be, Al, Fe, Cu, Ag, W, Pb, and U by magnetically 
analyzing the degraded proton energies of the 742-Mev proton beam 


of the Berkeley cyclotron after passage through the absorbing elements. 
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RANGE-ENERGY MEASUREMENTS FOR 742-Mev PROTONS 
John E. Ward 
Lawrence Radiation Laboratory 


University of California 
Berkeley, California 


Mies 4, 1959 


I. INTRODUCTION 


If a charged particle passes through a substance with sufficient. 
energy that the valence-sheill effects--i.e., atomic binding--can be 
meclected, the average rate of energy loss due to ionization onlyas 


expressed by the well-known Bethe equation, 


4 2 2 
_ GE _ 4re iz wo Cc go nies 4 ; (1) 





dx Zi I 


TLV 


mene ec and m refer to the electronic charge and mass, wezelogea. 
charge of the incident particle, NZ is the number of electrons per 
unit volume of the stopping element, and I is the 'mean excitation 
potential'' of the atoms in the stopping element. 

The development of this formula makes use of the Born approxi- 
mation, which requires the orbital velocity of the electrons to be small 
compared with the velocity of the incident particle. A proton witha 
kinetic energy greater than 2 Mev satisfies this condition. 

It can be seen that the stopping power depends only on the 
velocity and charge of the incident particle and not its mass, and the 
constant of primary importance as far as the stopping substance is 
concerned is its mean excitation potential. The mean excitation 
potential is defined as the product of all excitation levels in the sub- 
stance, Ao.» each weighted exponentially by its respective oscillator 
g 


strength, fio: 


fik 
I = HiT wo. 9 (2) 


where 


This formula would offer a means of calculation of I, but, except for 
the lightest of elements, a very impractical one. I is usually taken 
as an experimentally adjustable constant. 

Bloch has shown, using the statistical model of the atom, that 


in the Thomas-Fermi approximation one has 


I= KZ, (3) 


where K is an adjustable constant. : Using the more accurate Thomas- 
Fermi-Dirac model, Jensen found 
I = Ko ( 


The index 0 refers to isolated atoms. This function is plotted in 


L+k)z72/3) Z. (4) 


Fig. 1 as the solid line. 
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Fig. 1. Experimental mean excitation potentials and the 
statistical atomic model. O: I/Z of bound atoms meas- 
ured with 340-Mev protons and evaluated relative to 
I = 165 ev; @: Ip/Z of isolated atoms. A: i) Z nee 
ured with 10-20-Mev protons and corrected for inner 
shells. The corrections apparently underestimate the 
stopping -power deficiencies of inner shells at low 
energies. 
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Tl. EXPERIMENTAL 


Since at high energies the correctiong applied to the Bethe 
formula for valence-shell effects are negligible, the 742-Mev protons 
from the Berkeley cyclotron seem to be a natural choice for the incident 
particles. In addition, this energy range would serve as a check on the 
340-Mev Bakker and Segré data. z 

Figure 2 is a reproduction of a slide used by Werner Brandt at 
the Gatlinburg Meeting on Penetration (September, 1958), and compares 
the effective mean excitation potentials (without corrections), meas- 
ured by Burkig and MacKenzie at UCLA in 1957, ° with the fluctuations 
of the electron densities in the periodic system. One can notice slight 
but significant trends inthe I values, which are parallel to the electron 
density wherever the accuracy inI is high enough. The second max- 
imum should have been shown clearly by Mo(42) and Rh(45), but un- 
fortunately these were the poorest foils and hence gave the least accur- 
ate data inthe series. Lead showed an "anomalously" high value and 
needs an independent check. Thorium, onthe other hand, confirmed 
the expected trend. 

On the basis of this diagram and on the advice of Dr. Brandt 
the following metals were chosen as target elements: Be(4), Al(13), 
Fe(26), Cu(29), Ag(47), W(74), Au(79), Pb(82), and U(92). It was 
decided to use enough absorber to degrade the protons approximately 
15% in momentum or from 742 Mev to 570 Mev in kinetic energy. In 
each case an equivalent amount of absorber was to be used so that the 
results could be compared ona relative basis with aluminum as the 
standard. (See Table I) Adequate information exists on Cu, Al, Ag, 
Pb, and Be to determine the amount of absorber required. Specifically, 
the tables and graphs in UCRL-2426 (rev. )! were used. To determine 
the equivalent thickness of Fe, W, Au, and U a numerical integration 
of Bethe's formula was performed using values of I previously de- 
termined by Bakker and Segre (see Appendix A). Clearly these values 
are not precise. To provide for a means of interpolation in order to 
use the precise equivalent amounts of absorbers for the relative 


comparison, an additional length of aluminum absorber is to be added 
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Fig. 2. Mean excitation potential compared with electron 
density. (From Werner Brandt, Gatlinburg Meeting 
on Penetration, Sept. 1958. ) 
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to each absorber in order to provide an additional point. This absorber 
osham, 9 a fssal of aluminum, is to be placed first on the "upstream" 
and then "downstream" side of the various absorbers and the results 


averaged to give the effect of adding this additional absorber in the 


middle. The equivalent values of the aluminum shim are listed in 


Table II. 


Table J 








Equivalent absorber thicknesses 


a a, 


Aa ! 
mass stopping power 


Absorber Density Length 
(g/cm?) een) (cm) (in. ) 
Be(4) i235 85 45.94 18.09 
Al(13) 200 89 a2e9S ZS 
Fe(26 fies: 91 Pl 59 4.56 
Cu(29) S205 101 11.30 4.45 
Ag (47) 10.49 110 10.49 4,13 
W (74) eh oe iS 6.20 2.44 
Au(79) 18.88 118 6.24 2.46 
Pb(82) 11.34 esis! 1 Merges: 4.62 
Ut? 2) ro. 70 1Z5 6.68 2708 
Table II 
Equivalent value of Al shim 

Absorber Equiv. value* Thickness (g/cm*) 
Be (14) 1.024 9.216 
Al(13) 1,000 9.000 
Fe(26) 0.906 8.154 
Cu(29) 0.875 UsatD 
Ag (47) 0.789 fon 
W (74) 0.680 6.120 
Au(79) 0.676 6.084 
Pb(82) 0.660 5.940 
U(92) 0.630 5260.0 
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The experimental setup is shown in Fig. 3. The full energy of 
the 184-inch cyclotron is used, 742 Mev. The proton beam passes 
through a pair of jaws reducing it to 1/16 in. in the horizontal dimen- 
sion. The beam then passes through focusing and bending magnets 
and comes to a focal point 1 foot within the physics cave. This point 
is the physical center of the absorbers. The beam will then pass 
through a focusing magnet (three-element quadrupole), and a bending 
magnet in which the beam will be momentum-analyzed, and then fall 
upon x-ray film placed at the focal point of the quadrupole. 

The procedure will be to obtain an x-ray photograph of the full 
beam and then a photo of the beam with the absorber in place for each 
of the elements listed. In addition, photos will be taken with the 
aluminum shim as previously explained. 

The density of the developed x-ray film will be proportional to 
the number of protons striking it. By reading the film with a den- 
sitometer one can obtain a plot of the numbers of protons versus their 
positions on the film. In order to facilitate density measurements a 
duplicate of each of the exposures mentioned above will be made at 
precisely one-half the time for the original. Then, for example, the 
maximum intensity on the second exposure will be exactly one-half 
that of the maximum intensity on the first, and can be used as the 
reference for finding the points of half-maximum intensity on the first 
film. This procedure can then be extended ad infinitum throughout 
the intensity spectrum. The plot of numbers of protons versus film 
positions is expected to look like Fig. 4. In addition, by the use of 
wire orbits, a plot of energy versus position on film can be made This 
plot should look something like Fig. 5. The two graphs can then be 
combined to form a plot of number of protons versus energy. An ex- 
ample of this type of plot is shown in Fig. 6. In general it has the same 
shape as Fig. 4. This plot (Fig. 6) will be the tool for describing the 


experimental results. 
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Pig. 2. Experimental setup, showing placement of magnets. 
Net to scale. 
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Fig. 4. Expected general contour of plot of the number of 
protons, i.e., the film density, versus position on an 
x-ray film after undergoing magnetic analysis. Two 
curves are shown: that of the full-energy proton beam 
of the 184-inch cyclotron, 742 Mev; and that of the 
same beam after approximately 15% momentum 
degradation. 
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Pig Owe oxpeeree oener al contour of the results of a wire- 
orbit analysis of an analyzing magnet. 
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Fig. 6, Expected results of the combination of data ob- 
tained from figures 4 and 5 showing the energy spectrum 
of the full-energy proton beam, 742 Mev, and the same 
beam after approximately 15% momentum degradation. 
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Il, DATA ANALYSIS 


For determination of the mean energy loss through each 
absorber the center of gravity of the two curves of Fig. 6 must be 
found. This can be done easily for the first curve, representing the 
full beam of protons undegraded. However, because of close collisions 
of protons and electrons the second curve has, theoretically, an 
infinite tail. Obtaining this infinite tail is not practicable experimentally, 
so that it must be supplied by a theoretical calculation. The observed 
ranges of individual particles from any monoenergetic source will show 
a substantially normal distribution about the mean range. Because the 
very hard collisions are few in number, the actual distribution is some- 
what asymmetric, with a long tail in the direction of the short ranges, 
and a mean range which is slightly shorter than the modal range. : For 
the approximation of this tail, the collisions suffered by the protons 
with the electrons of the absorbers will be divided into two classes, 
"near" and "far. '' The so-called near collisions are those responsible 
for the occasional large energy losses and the far collisions are those 
contributing only to the mean energy loss, but introducing negligible 
energy-loss fluctuation. In other words, each proton will suffer 
approximately the same minimum energy loss and may in addition 
suffer a near collision (or collisions) resulting in a much greater energy 
loss. Each absorber will then be divided into a number of equal seg- 
ments in which the probability of a near collision is known, the energy 
loss will be computed for this one segment, and the resulting spectrum 
will be "folded" into itself as many times as there are segments. It 
is hoped that this spectrum will fit the experimental curve and provide 
the needed tail. This procedure seems valid in that the so-called near 
collisions are those in which the binding energies of the electrons are 
insignificant, and this is the assumption made originally by neglecting 
Peeweorrections to Bethe's formula. 

By this means of completing the graph for the degraded protons, 
the center of gravity of the curve can be found and the mean energy loss 
determined. This procedure would be done for each of the elements. 


By interpolation between the two curves, obtained by use of the 
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aluminum shim, for each of the elements the thickness of absorber 
equivalent to aluminum will be determined. The relative mass stopping 
power of each element will be calculated by dividing the amount of 
aluminum (in aera by the amount of the element (in ese) needed 

to produce the same mean energy loss. 


In addition, the stopping power per electron is, from Bethe's 





formula, 
dE 
relax 7 aeeone Dray 2 Shea at 
SNS Sree 2 I eae ; 
mv 


The relative stopping power per electron with respect to aluminum is 


Serie (eo oe Sy riglecte ac, 


= 


- gnI gn 2mv~ - £n(1-8%)-B*] 


ag tl 


This quantity can be easily calculated by multiplying the relative mass 


stopping power by - 
Ax 


a ee 
A Z 


Af 
Thus the means for calculating I is apparent. In this calculation the 
value of I = 150 ev for aluminum is used. 

From the values of I calculated, a curve of 1 Z versus Z will 
be plotted. This curve will in general look like Fig. 1. The value of 
the Bloch "constant" can then be compared with previously obtained 
valWes. 

In order to determine the length of absorber traversed by the 
protons a correction for multiple scattering must be made. In the 
usual approximate correction, the observed range is expressed as a 


projection of the mean range: ae 


eee 
ers 2 #;c08 8; > y 2,01 = 5 7) 
1 1 


=e ee 7 e.* , 
mean 2 1 1 
al 


By | 


where f. is the distance between the ith and the (i-1)st small- 
angle collision and 0. is the direction with respect to the beam axis 
yy 


after the ith collision. Therefore we have 
1 Robs 5 
mean ‘obs “~~ B ( os (x) ) oS 
0 


where (6, (x)) is the mean square deviation in angle from the 


normal direction and x is the projected distance in the absorber. 


Also, 
Sema Vg) = BME 


where t is the absorber thickness measured in radiation lengths. = 


Using the einer 


0 
(0° (x) 4 BNE <2 (Gov in dee 
x 


Bichel, Mozley, and Aron domme 


poiachd ame G12 Ro1./R 


mean mean 


Eo 1,2 


Gm m7, Uf “0 p° dE dE 
0 0 B/Z B<(B/Z) 
2D E 2 
(M ) 
0° f az > 4d 


Where G= 2£n181 z-V/3. B is the stopping number of the atom, 


2 
B=Z c eS. 6 | | 
ep) 


and M4 is the rest mass of the electron. The ratio AR/R is a slowly 


varying function of the energy, amounting to several percent for high Z. 


oe li 


The corrections found by these experimenters appear in Table III. In 
each case it is clear that as Z increases the correction for multiple 


scattering increases, andas E increases this correction decreases. 





Element E Range 1s 
(Mev) Path length 
Be 17.340 0.9790 
9.588 Or 733 
Al 17.836 0.9962 
14.971 0.9960 
LisZ0 0.9956 
6.150 009925 
Cu 17.893 0.9900 
SP SIETS: 0.9872 
Ag 17.923 0.9818 
10.022 OB STOaS: 
Au We Aare) 0.9626 
9.698 0.9475 


“correction for multiple scattering 
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APPENDIX 


A. Determination of Absorber Thicknesses 
l. Fe, W, Au, and U 
Since no usuable data exist on these elements, a numerical 
integration of Bethe's equation was performed in order to determine 
the absorber thickness. The total desired energy drop was divided 
into twelve equal steps such that the change in B was constant, 


0.004. The energy increment was then calculated from 


Z 
E =MpyYec ; 


where 





Therefore we have 
= Za 2 _ 3 
AE = mM 9¢ By AB = 3.752 By (Mev). 


The integration on Bethe's equation was performed by de- 








termining © 
is a. <cppZ mp a 
ax/" 25 |) eee ee 
pA B G 
where 
4,2 4 
B= 4 tre = No  _ POT g cm atoms 
mM ¢ See g-atom 
ae wo 


The values determined appear in Tables IV and V. Graphs of the re- 
sultant variation in dE/dx are shown in Figs. 7 through 10. 
2. Be, Al, Cu, Ag, and Pb 

Sufficient data on these elements exist in UCRL-2426 (rev). £ 
a. oul ey. 


The absorber thicknesses to be used were given in Table I. 
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bag. (2 ‘Calculated range=energy curve for He. 
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Fig. 9. Calculated range-energy curve for Au. 
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B. Bending-Magnet Analysis 


1. Angle of Bending 
If the H£ of the magnet and the Hp of the particle are known, 


the following is the derivation for the angle of bending (See Fig. 11): 


— 
ae. (2p 


3 


_ ee patericlle pane, r 
= ~~ (AD) 
magnet mag. 
sin g = mae hag. : 
etene Dart. 
(Hf) 
ae sin J = then §@ = Eo 
(FHP) part. 


This approximate formula works well for 09 less than 1 radian. 
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iia 
fe) 


Fig. li. Path of a particle in a magnetic field. 


a. 


é. Vertical Focusing in an Analyzing Magnet 





Fig. 12. Side and top views of typical bending magnet, showing 
field configurations and entrance and exit angles. 


Vv 








In Fig. 12, let 0, = 0, = 9 and tan@,=-—~ . From Maxwell's 
l 2 Z 1 Vv 
uy 
equations we have 
VXH =0, 
0H 0H 
a I ef) 
OX oy ‘ 
dH 0H 
en 2 =U 
dy OZ ; 
0H dH (1) 
x Z 
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I shall let He = Omethnen 





aH oH aH. 
a le 2) 


and substituting Eqs. (2) into (1), one has 























dH dH dH dH dH 6) lal 
== = O a Ms ee z 
ax ox ; dy oz y OZ 
Using Gaussian units, 
F = (e/c VXH 
yoo / i 
= e/c(v H er) = -e/c v, H,. 
oH, 
Assuming ee is constant leads to 
OH. oH wee 
dy a dy Y" 2 
0 
Therefore 
oe 
F =-e/cv_H_ =-e/c vy oe OY 


oe 


The momentum imparted to the particle inthe y direction is 


q ee dH 
FP dt2@ fe =e) 22a 
y ye € a dz 


= a eS 
- tan Oy dizer = tand,y ATi 


Ap 


dz yi 


Defining 5 = Ap,/P: where p =eHp/c, one has 


(e/c)tan 0, yOH, 


OS “e/cjHp = tan 0, 


ol 
Md 


In the pole gap H = Be and ek is the increase of a from a place 
of zero field to the edge of the pole tips, where Be becomes uniform 


and equalto H. Therefore 


and 
e 3 
p (3) 


= 30) — 


I have assumed oy small until reaching the pole tips, so that 


ve 1s constant. From Fig. 13 we have 


y = S oP = Sa. 


Pie. lo. Vertical bending at pole tips due to field variations. 


Therefore 
= a ae yee 
S6=at+B = = + é 
1 0 
or 
i) 1 § 1 
a ee ee : 
S69 S. y ate 
from iq. (3), 
Ab = 6 2 tan 0) 
f y p 


This means that f °*1s independent of the height at which the particle 
hits the entrance of the magnet. This independence contributes the 
focusing action in the vertical plane. The principal planes are 
assumed to be the ends of the magnet. In terms of the total angle of 
bending, 


iI Z 
= = (4) 


a3 e 


Change in Vertical Focal Length with a Change in Momentum 
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4. Dispersion in an Analyzing Magnet 

In order to determine the dispersion particles of various 
momenta passing through an analyzing magnet, it is necessary to find 
distance d in Fig. 14. Figure 14 shows two particle tracks, with 
point (h, k) being the center of curvature of the lesser momenta 


particle. From the geometry we have 


h=- L/2+ p! sin(@/ 2), 
k = p' cos(6/ 24 
The equations for the circles s and s' are 
(x - 0)* + (y - pcos(0/d)* = p*, | 


(x + = - p' sin(@/2)* + (y-p! cos@/2)° = p' or 


To find d it is then necessary to find the intersection of the circle 
s' and the line x = L/2. Note that we have 
J Ze 
L/2 cos(@/2) = ea a [4 


in(0/2)= ——, 
sin(8/2)= = 5 


Therefore the equation for the circle s' is 


[x+(p-p' ) sin@/2)] 4 (y -p' cos(6/)° =p! é 


— 2x(p-p') sin(@/2)+ (p-p' yo sin“(0/2)+ a ~ 2yp' cos(@/2)+ (p' \-cos 10/ 2m p' 


a + 2x(p-p') sin(0/2)+ mean 2) Z2pp' sin“(6/2)+ e - 2yp'cos@/2)= 0. 


Now at x = L/2 (i.e., the intersection of circle s' and line x = L/2), 


we have 


2 
= + L(p-p' ) sin(@/2)+ oe sin“0/2)- 2pp' sin“(0/2\+ yo - 2yp' cos(6/2) = 0. 


This reduces to 


= p' cos(6/2) - J oy * cos (0/2) - fie Oo) : 


= 040 
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Fig. 14. Dispersion in an analyzing magnet of charged 
particles of different momenta. 
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Since p=L/2 sin , 


we have 
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C. Wire Orbits 
l. The Wire-Orbit Method P 
This is an experimental method for determining the path ofa 
charged particle in a magnetic field. It is based on the analogy be- 
tween the shape assumed by a current-carrying wire in a magnetic 
field and the path of a moving charged particle in the same field. 
The motion of a charged particle in a magnetic field is described 


Dyethe ecauation 





Vs 
mv _ > —> 
R = @y <5 
with the solution 
BR = 72 (1) 


The forces on the wire in Fig. 15 consist of the tension, T, 


and the force due to the magnetic field, 


— 


Feidé xB. 
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Fig. 15. Forces involved in the wire-orbit method. 


Equating the forces, 
2 sind = 1 dae, 


and noting 


dz/2 


R as df > QO, 


Sit GQ = 








eo). 


one finds 


Eee (2) 


i 
Therefore it 1s clear that if the ratio of tm sion to current in the wire 
1s the same as the ratio of momentum to charge of the particle, the 
wire will assume the path of the charged particle. 

The above derivation is for a constant B, but it is possible to 
divide the wire into many small segments where the immediate field 
1s considered constant, and, in the limit, the wire orbit and particle 


path will be the same for 


— 
_— 


al 
a 


0 fro 


(3) 


In practice the wire should be very fine and flexible so that the 


force of gravity may be neglected. In workable units, 


p(Mev/c) = aS , 


where M indicates the weight attached to the wire to provide the 
tension. 


2. Example 


Tie 
point Wires 





Fig. 16. Experimental setup in wire-orbit method. 
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A typical experimental setup for using the wire-orbit method 
is shown in Fig. 16. The two pieces of wire passing through the 
magnet are in series so that they represent two possible "rays" of 
particles eminating from the tie point, A. Actually the tie point must 
be located slightly ahead of the planned point of particle emission, 
since the wires occupy physical space and cannot be made to emanate 
from a point. The wires then pass through the magnet, past a frame 
for measuring their separation, and then to a point just short of a 
pulley. A short piece of insulated wire is tied to the current-carrying 
wire and then passed around the pulley. This is done so that the 
Cumments injeach wire have the same direction sense.  Usimpo ities 
point as the object and the intersection of the wires, as determined 
by a projection of the measured separation at the frame and diameter 
of the pulley (as in Fig. 17) as the image, one solves the focusing 
problem by simple geometry. In practice the current in the magnets 
is adjusted until the separation, s, at the frame is correct to provide 


the desired image point. 


wires 





separation at frame 
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pulley diameter 


Fig. 1?. Geometric determination of focus. 


aS Ca 


FO, 
ie 


ie. 
(So 


e302 


REFERENCES 


M. S. Livingston and H. A. Bethe, Revs. Modern Phys. 9, 285 
ge). 

W. Brandt, Health Physics 1, 11-20 (1958). 

Meecloch,) Z. Physik 81, 363-376 (1933), 

H. Jensen, Z. Physik 106, 620-632 (1937). 

C. J. Bakker and E. Segre, Phys. Rev. 81, 489-492 (1951). 

V. C. Burkig and K. R. MacKenzie, Phys. Rev. 106, 848 (1957). 
Jo oH. Atkinson, Jr. and B. H. Willis, High-Energy Particle Data, 
Vol. II, UCRL-2426 (rev), 1957. 

D. C. Sachs and R. J. Richardson, Phys. Rev. 83, 834 (1951). 
R. R. Wilson, Phys. Rev. 60, 749 (1941). 

R. Mauer and E. Segré, Phys. Rev. 84, 191 (1951). 

B. Rossi, High-Energy Particles, (Prentice-Hall Publishing Co., 
New York, 1952). 

B. Rossi and K. Greisen, Revs. Modern Phys. ish 240 (1941). 
Bichel, Mozley, and Aron, Phys. Rev. 105, 1788-1795 (1957). 





This report was prepared as an account of Government 
sponsored work. Neither the United States, nor the Com- 


mission, nor any person acting on behalf of the Commission: 


ANE 


Makes any warranty or representation, expressed or 
implied, with respect to the accuracy, completeness, 
or usefulness of the information contained in this 
report, or that the use of any information, appa- 
ratus, method, or process disclosed in this report 
may not infringe privately owned rights; or 


Assumes any liabilities with respect to the use of, 
or for damages resu]ting from the use of any infor- 


mation, apparatus, method, or process disclosed in 
this report. 


As used in the above, "person acting on behalf of the 
Commission" includes any employee or contractor of the Com- 


nmassion, 


or employee of such contractor, to the extent that 


such employee or contractor of the Commission, or employee 
of such contractor prepares, disseminates, or provides access 


to, 


any information pursuant to his employment or contract 


with the Commission, or his employment with such contractor. 























